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Abstract 

We propose schemes to create cluster states and W states by many superconducting-quantum-interference-device (SQUID) 
qubits in cavities under the influence of the cavity decay. Our schemes do not require auxihary qubits, and the excited levels 
are only virtually coupled throughout the scheme, which could much reduce the experimental challenge. We consider the cavity 
decay in our model and analytically demonstrate its detrimental influence on the prepared entangled states. 
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Entanglement is an essential resource for testing quantum nonlocality and also for quantum information process- 
ing. Although many proposals have been put forward, the achievement of multi-partite entangled states are still 
challenging experimentally due to decoherence as well as the limitation of current techniques. This work focuses on 
a combinatory system with superconducting devices embedded in a cavity. Among a variety of qubit candidates, 
superconducting devices such as Josephson junction circuits 1, 2J, Josephson junctions [MBQIj Cooper pair boxes 
0,01, and superconducting quantum interference devices (SQUIDs) [1, S Hfl > have drawn particular interests 

due to their potential for integrated devices. On the other hand, cavities are currently considered as excellent devices 
to achieve niultiqubit entanglement for atoms, ions, quantum dots, and other charge qubits 13, 14, 15, 16]. 

Our goal is to present ways for creating multi-partite entanglement of many rf SQUID qubits in cluster states 
and in W states in cavity QED. The cluster states is the key ingredient in a measurement-based quantum 
computing (QC), i.e., one-way QC [3i with which we can carry out expected quantum gates by some singlc-qubit 
operations and detections. The general definition for the cluster state |0{fe})p is from the set of eigenvalue equations 

K^°''> |0{fc})(^ = (— 1)'^" ' with K'^"-^ ~ cri"'' (g) cri''\ and nghb(a) specifying the sites of all the neighbors of 

h£nghh(a) 

the site a and Ha G {0, 1}. It has been shown that the cluster states can be created by photons in linear optic system 
and by atoms going through cavities Ii9, 20]. While for doing a one-way QC, the flying qubits, such as the moving 
atoms and the photons, are not good candidates in view of an accurate manipulation on QC. In contrast, SQUID 
qubits embedded in cavities are always static, which would be good for one way QC if we could prepare them into 
cluster states. W states, whose general form is Wn — {l/\/n) \n — 1, 1) with \n — 1, 1) being all the totally symmetric 
states involving n-1 zeros and 1 one, are famous for their robustness to local measurement, even under qubit loss. 
This makes W states useful for, e.g., quantum communication based on many nodes of a network. A lot of schemes 
for W state preparation by atoms going through cavities have been proposed, e.g., a very recent scheme for rapidly 
creating multi-atom W states under the cavity decay [IJ. Compared with the moving atoms, the SQUID qubits, 
without movement in the cavities, are obviously more suitable for this job. 

One of the favorable features of our schemes is that we only virtually couple the excited level. So our implementation 
subspace is only spanned by the two lowest levels of the SQUIDs, which makes our scheme simpler than previous 
work llCtllll . and very robust to decoherence due to spontaneous emission from the excited level. Moreover, auxiliary 
qubits [20j are not required in our scheme, which improves the efficiency of our implementation. Furthermore, we 
consider in our treatment some imperfect factors, such as the cavity decay and various coupling strength of different 
SQUID qubits to the cavity mode and the external microwave. 

For each SQUID, two lowest flux states and an excited state are employed in our model, as shown in Fig. 1 
where quantum information is encoded in ]0) and ]1). The spacing between two neighboring SQUIDs is assumed 
to be much larger than the size of each SQUID ring (on the order of 10-100 /im) so that the interaction between 
any two SQUIDs is negligible. The SQUIDs are radiated by the cavity field and by a classical microwave pulse 
which are tuned respectively to ujc = (wo2 — S) and cj^^j = (wi2 — <5), with wo2 and Wi2 resonant frequencies between 
levels ]0) and ]2) , and ]1) and ]2) , respectively. S is the detuning. The Hamiltonian of the SQUID can be written 

in an usual form 0,112, i.e., Hsj — -^f- + Ej cos ^27r^^^ , where the conjugate variables and Qj 

are, respectively, the magnetic flux threading the ring and the total charge on the capacitor of the jth-SQUID, with 
the commutation relation [^j,Qj] = ih. is the static external flux applied to the ring of the jth-SQUID, and 
Ej = Icj^o/'^Tr is the maximum Josephson coupling energy of the jth-SQUID with Icj the critical current of the 
junction, and $0 = h/2e the flux quantum. Under the action of a single-mode cavity field He = uJc{a^a + i), with 
H = 1 assumed, and and a the creation and annihilation operators of the cavity mode, we have the interaction 
Hamiltonian Hj, Hj = X^jLi ^T'(^J ^ ^xj)i^cj + ^fj.wj), where $cj" and ^f^wj are the magnetic flux threading the 
ring of the jth-SQUID generated by the magnetic component 'Bj(~r,t) of the cavity field and the microwave pulse, 
respectively. $y — 'Bj{~r,t) ■ dS j with i = c, iiw corresponding to the cavity field and the microwave pulse, 
respectively, and j — 1,2,.. . , N regarding different qubits. 

Following the steps in [l^], we get to the Hamiltonian of the system in the interaction picture = 
'Ef=i gjci^e-'^^ \0).. (2] + .gj-ae'^* 12)^.^. (0] + %e-'"^* Jl)^.^. (2] + Qjc'^* \2).. (1] .We can adiabatically eliminate the ex- 
cited level ]2) by a similar way to in 23] and reach an effective Hamiltonian, Hj = J2f=i ("^^ (1] + a Wjj i^lj j 

where Xj = and Max {g^, r2j}<C 5j. Considering the influence from the cavity decay, we obtain following 

Hamiltonian, 

= ^ A, (at 10)^.^. (1| + a 11)^.^. (01) - z|ata, (1) 



2 



where k is the cavity decay rate, and we have assumed that k is much smaller than so that no quantum 
jump due to the cavity decay actually occurs during the time evolution under our consideration. To obtain clus- 
ter states, we assume that the N SQUIDs are prepared initially in a product state \^po) = (^f=i^ \+)j |0)^ where 
1+)^ is the eigenstate of al with eigenvalue 1. The cavity field is prepared in (|0)^ + i , which is from 

a resonant interaction of an ancilla qubit initially in the state (|0) — |1)) /^/2 with a vacuum cavity. Before we 
get started, we assume the SQUIDs to be decoupled from the cavity field, and the microwave to remain off. By 
adjusting the level structure of the SQUID 1 and turning on the classical microwave to meet the condition in 
Fig. 1, we reach Eq. (1) with N=l. A straightforward solution of Eq. (1) could yield |0)^|0)^ |0)c , 
|0)Jl)i - e-^{[4§^sin(Gii)+ cos(Git)]|0)Jl)i - *ALsin(Gii)|l)jO)J, |l)jO)i ^ sin (dt) + 

cos(Git)]|l)jO)i -z^sin(Gi<)|0)Jl)J, ^ e-^\l)Jl),, where Gi = ^A^ - f|. If we choose the 

evolution time of the system to be ti — [arctan(— + 7r]/Gi, the state evolution of the system is given by 

lV'>i = 5[|0)c(|0)i+e-'^^sin(Giii)|l)i)+z|l),al (e^^ ^ sin (Giii) |0)i + e"^ |l)i)] ® J] l+>, |0)Ar, where 

al = (1| — \0)ii (0| , and the state is not normalized if k is not zero. The cluster states prepared below will be 
written also following this convention. 

Our next step is to adjust the level structure of the SQUID 1 back to the previous situation with decoupling to the 
radiation field, but adjusting the level structure of the SQUID 2 to do the same job as done on the SQUID 1. So we have 
= i7l{|0)e [10)2 +aie-^ ^ sin (02^) \yi)]+ 1 ^^[e-^^ sin (Gaia) lO)^ \x^)+ale-'^ \y{)]}® 

N-l 

n where \xi) = |0>i+ e-^^ sin (Gi^i) , = e"^ ^ sin (Giti) 10)^ + e"^ and G2 - 

i=3 

\/'^2 ~ ft effective coupling regarding the SQUID 2, including the cavity decay. Step by step we may carry out 

above operation on the rest SQUIDs individually. After the (N-l)-th SQUID is performed by the same operation as 
on the first one, we have 

I^)a^-i ^ ^(lo)c I^A'-i) + * Ii>c^f \yN-i)) ® |0)^ , (2) 

where \xn~i) and \yN-i) are two different entangled states of the first (N-1) SQUIDs, and they follow re- 
cursive relations below \xn-i) = |0)^_ Ja;Ar-2) + erf ^-^ gin (GAr_iiAr_i) |1)^_ Jy7V-2) , \yN-i) = 
e-^^T^^^sin(GAr_itAr_i) |0)^_;^ |a;7v-2> + crf~^e"^^^^ |l)Ar-i I2/JV-2) • For the Nth SQUID we stop the state 
evolution at ijv = cJ^[8'i^ctan(— ^^^) -|- tt], with Gn — \J ~ f|: then can obtain 

l'^)^ = ^ (lO)iv I^A^-i) + ^^'"^ {GNtN) |1>^ <7^-' \yN-i?j ® |0), . (3) 

To check above state, we first assume k = 0, i.e., the ideal condition. In this case, Eq. (3) reduces to the standard 
cluster state [l7j . 

i'f)A^-^<^i(|o). + |i).-rO- (4) 

To evaluate our prepared cluster state under the cavity decay, we have calculated the fidelity F and the success 
probability P, as shown in Fig. 2 where we considered the qubit numbers 2, 3, and 4, respectively. It is evident that 
both F and P are going down with the qubit number and k. It is understandable that the preparation of a cluster state 
with bigger size need longer time and is thereby affected more by the cavity decay. Thus in the presence of the cavity 
decay, the size of the cluster states are limited. For example, if we want F > 0.95, in the case of Ai = A2 = ... = A, the 
maximal qubit number of the cluster state is 32, corresponding to F = 0.951. For a cluster state with more qubit, we 
could created it step by step by connecting the few-qubit cluster states by quantum phase gates on the end qubits 
of the neighboring SQUID chains Moreover, to keep the system as stable as possible, we should not change the 
level spacing of each SQUID drastically. In this sense, we estimate the appearance of imperfection due to the time 
delay of the level spacing adjustment, i.e., the undesired phases regarding the component states of the superposition. 
This estimate is plotted also in Fig. 2. 

From now on, we turn to the generation of W states using SQUIDs. In contrast to the cluster state preparation by 
individually manipulating the SQUIDS, we have to simultaneously couple all the SQUIDs in this case to the cavity 
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mode and to the microwave. Provided that the SQUIDs and the cavity mode are initiaUy in the state Y[ |l)i |0)i ^^^d 
the vacuum state \0)^, respectively, we have the time evolution of the system by straightforwardly solving Eq. (1), 

m)) = n |0), |0), {1 + ^hl + e-^(cos ^ + ^ sin ^)]} 



En \ T-T ^l^k r ^ '^^ / B^ ^ ■ B^ \i 

\W n 10), 10). ^[-1 + e-- (cos - + - sm -)] 

^ iiXi . Bt - 



- njO),|l),-gisin— e-T. (5) 

N 

where = J^^p and B ^ \/l6A2~^. In order to obtain W states of (N-1) SQUIDs, we have the conditions 
t = and = A'^e^' to be satisfied. Then the system evolves to 



N 



M^A^-i = e-^^^|l,) n |0),, (6) 

N 

where A'^ = ^ X^- This is a, W state of (N-1) SQUIDs with arbitrary coefficients. By choosing the same coupling 
for the rest (N-1) SQUIDs except the first one, we may produce a standard W state 

t 1 A ^ 
W^^-i = e-^^==^|U.) n |0),, (7) 

with the success probability P = e^^. We plot in Fig. 3 the success probability P versus cavity decay k. From 
above condition for the evolving time, we have Ai = A'e^'^^^-^ > A', which means that for many-qubit case, the 
coupling strength regarding the first qubit is much larger than those of others. This also implies that, the coupling 
strength regarding other SQUIDs would be quite small if Ai is not big. Experimentally, the difference of Ai from other 
coupling strength could be made by reducing the detuning regarding the first SQUID. Alternatively, we may put the 
first SQUID on the antinode of the cavity standing wave (i.e., the maximum coupling), but others on the positions 
deviated from corresponding antinodes. 

We address below the experimental feasibility of our schemes. The implementation time of our scheme should be 
much shorter than the cavity decay time — Q/2tti/c, where Q is the quality factor of the cavity (Q= 10^ — 10^ 
has been achieved experimentally [2^, and for superconducting qubits [ISl), and is the cavity field frequency. 
The coupling constants of the SQUID to the cavity field and to the classical microwave available at present are 
g ~ 1.8 X 10^s~^, and ~ 8.5 x 10^s~^, respectively [23|. So we may have an effective coupling A ^ 10 MHz if 
assuming S ^ 1.5 GHz. In this case, the interaction time r is on the order of 10~^ sec, much shorter than Ax 10~^ 

(with Q — 10^ and fc = AOGHz). This implies that our proposed states can be of the fidelity higher than 0.99 due 
to k/A 4 X 10-2. 

Although optical photons in cluster states have been experimentally demonstrated some important features for 
one-way QC, we prefer fixed qubits, like the SQUIDs embedded in the cavity, for storing and processing quantum 
information. This statement is also applicable to the W state case for quantum information processing. Moreover, 
due to no necessity of auxiliary qubits, our scheme could much reduce the experimental challenge for cluster state 
preparation. For example, generation of a two-qubit cluster state requires at least five steps in 19] due to the auxiliary 
states involved, while only two steps are needed in our scheme. Furthermore, we only virtually couple the excited 
levels throughout our scheme, which, compared to [lOllllj . could both simplify the implementation steps and reduce 
the possibility of decoherence due to spontaneous emission. We have given the analytical expressions of the prepared 
cluster state and W state under the cavity decay, and numerically assessed the related fidelities. 

Before ending our discussion, we give few comments on our wavefunction treatment for the cavity decay. The 
quantum-jump approach has been widely employed in solving dissipative dynamics for quantum systems, which 
lies in the time evolution governed by a non-Hermitian operator and interrupted by instantaneous jumps by the 
detection of a photon. Normally, a solution by quantum-jump approach has to be resorted to numerical calculations. 
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In our case, however, to prepare high-fidelity cluster states and W states, we prefer our implementation time much 
shorter than the cavity decay time. Therefore, we may only focus on the time evolution of the system governed 
by the non-Hermitian Hamiltonian Eq. (1) before the leakage of photons occurs. In this sense, our wavefunction 
method would be advantageous over the quantum-jump approach in the availability of analytical solution clearly 
demonstrating the detrimental influence from the cavity decay on the prepared states. So our solutions might be 
helpful in experiments for estimating the infidelity and correcting the error, particularly for implementation with 
cavities of comparatively low quality. This also implies that our numerical results in the last two figures are valid only 
within the regime of k/X <C 1 {e.g., k/X ^ 0.1). 

In summary, we have proposed potential schemes for creating cluster states and W states of many SQUIDs. Fast 
adjustments of the level spacings of individual SQUIDs are needed in the generation of the cluster states. As no 
auxiliary qubits or flying qubits involved, our scheme gives good candidates for one-way QC. In the generation of the 
W states, all the SQUIDs are coupled to the radiation fields simultaneously, which results in that the cavity decay 
only affects the prepared states globally, instead of on the internal structure of the W state. Throughout our schemes, 
the excited level is only virtually coupled, and we have specifically studied the detrimental influence from the cavity 
decay. Our analytical results have shown clearly that the cluster states and the W states may be generated with 
high-fidelity only in the case of tiny cavity decay rate. 

This work is supported by National Natural Science Foundation of China under Grants No. 10474118 and No. 
60490280, by Hubei Provincial Funding for Distinguished Young Scholars, and by the National Fundamental Research 
Program of China under Grants No. 2005CB724502. 
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Captions of the figures 

Fig. 1. The level diagram of a SQUID with three lowest levels |0) , |1) , and |2), where ujc and oj^u, are frequecies of the cavity 

and the microwave, respectively, gj = —^^{hujc/2^o) ■ {0| $ \2). Bj(T*) • dSj and Qj = -j^ {1| $ |2)^. Bj(T^) ■ dSj 

are, respectively, the coupling constants of the jth-SQUID to the cavity field and the microwave. 5 is the detuning to the 
excited level |2) by the radiation fields. 

Fig. 2. The fidelity and the success probability versus the cavity decay in the generation of a cluster state, where the solid, 
dashed and dash-dot lines represent the consideration with the time delay Ta = 0, O.Olr, 0.05r for adjusting the level 
spacing, respectively, with r the desired time in the perfect case. The curves from the top to bottom correspond to N= 2, 
3 and 4. For simplicity, we have assumed the same coupling strength for each SQUID qubits to the radiation fields, i.e., 
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\i = \ (i=l,2, 3, ....). 

Fig. 3. The success probability versus the cavity decay in the generation of a W state, where we assume ~ 1.8 x 10* s 
Qi ~ 8.5 X lO'^s-^ and <5 ~ 1.5 x 10^5"^ 
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